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Top-antitop quark pairs belong to the most abundantly produced and precisely measurable heavy- 
particle signatures at hadron colliders and allow for crucial tests of the Standard Model and new- 
physics searches. Here we report on the calculation of the next-to-leading order (NLO) QCD cor- 
rections to hadronic W"''W~bb production, which provides a complete NLO description of the 
production of top-antitop pairs and their subsequent decay into W bosons and bottom quarks, in- 
cluding interferences, off-shell effects, and non-resonant backgrounds. Numerical predictions for the 
Tevatron and the LIfC are presented. 

PACS numbers: 12.38.Bx, 12.38.Cy, 13.85.-t, 14.65.Ha 



The top quark is the heaviest of all known elementary 
particles and is expected to play a key role in any theory 
of the flavour sector of elementary particles. Its precise 
investigation is, thus, of great importance at the current 
hadron colliders Tevatron and LHC, where top quarks 
are mostly produced via top-antitop (tt) pairs. 

The first step towards precise theoretical predictions 
for tt production at hadron colliders was made already 
about 20 years ago with the calculation of QCD cor- 
rections at next-to- leading-order (NLO) [ll-13|- Later 
also electroweak radiative corrections were calculated 
0], and recently important progress has been achieved 
both in the resummation of logarithmically enhanced 
terms and towards the inclusion of QCD correc- 

tions at next-to-next-to- leading-order (l^ - [2^ . 

The above-mentioned predictions are based on the ap- 
proximation of stable (on-shell) top quarks, i.e. the top- 
quark decays, which proceed into pairs of W bosons and 
b quarks in the Standard Model, were ignored. Re- 
cently also studies [ii-El] at the NLO QCD level have 
been presented that include the top-quark decays via a 
spin-correlated narrow-width approximation, i.e. the top 
quarks are still on shell. In this letter we present first 
results^ at NLO QCD on the further generalization that 
the intermediate top quarks can be off their mass shell, 
i.e. we consider the process of W'''W~'bb production, in- 
cluding leptonic W-boson decays. 

The reaction pp — > W+W^bb -I- X represents one of 
the few remaining 2 — 4 LHC background processes on 
the Les Houches wishlist [28[. While various such 2 — > 4 
NLO QCD calculations have been performed in the re- 
cent years (see e.g. Ref. [ll] for a review), W"'"W~bb 



production involves the treatment of resonant particles 
for the first time in a hadron-coUider environment on 
that level of complexity. The two resonances can be con- 
sistently treated in the complex-mass scheme that was 
introduced at the NLO level in the context of the cal- 
culation of the electroweak corrections to the processes 



e+e — )■ WW — )■ 4fermions 
full NLO calculation for a 2 



^ Similar results on WWbb production have recently 
been shown by the HELAC-OPP collaboration at the 
meeting of the HEPTOOLS network at Granada (see 
|http://indico.cern.ch/conferenceDisplay.py?confId=91923| . 



30|, which was the first 
4 particle process. 

At leading order (LO), hadronic W+W^bb production 
proceeds via partonic channels with quark-antiquark (qq) 
and gluon-gluon (gg) initial states. A few representative 
diagrams are depicted in Fig. [T] In addition to doubly- 
resonant (DR) diagrams, where the W"'"W~bb final state 
results from the decay of a tt pair, our calculation also 
includes singly-resonant and non-resonant contributions. 
As is well known, the bulk of the inclusive W"'"W~bb 
cross section is efficiently reproduced by the widely used 
narrow-width approximation, which incorporates all DR 
effects in the limit of vanishing top-quark width. Ft — >■ 0. 
By including all off-shell effects from doubly-, singly-, 
and non-resonant diagrams, our calculation consistently 
describes all contributions that are suppressed by one or 
more powers of Ft /mt • These extra terms are mandatory 
in order to achieve percent-level precision in the (inclu- 
sive and differential) description of tt production, and for 
a reliable simulation of off-shell W+W~bb final states. 
To describe top-quark decays in a realistic way we also 
include the leptonic W-boson decays W"*" — >■ and 
W~ — >■ t'^/J." in a spin-correlated narrow-width approxi- 
mation. 

In the following we briefly sketch the calculation of the 
virtual and real corrections. A more detailed description 
will be presented elsewhere. In order to prove the cor- 
rectness of our results we have evaluated each ingredient 
twice and independently. The treatment of the virtual 
QCD corrections to qq/gg W"'"W~bb is based on di- 
agrammatic representations of the one-loop arnplitudes 
and numerical reduction of tensor integrals 3l|, |32| . The 
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agator, and the top mass counterterm S^t is related to 



FIG. 1: Representative LO diagrams of doubly-resonant (up- 
per line), singly-resonant (first diagram in lower line), and 
non-resonant type (last two diagrams in lower line). 
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FIG. 2: Hexagon diagrams in gg/gg — > W^W bb. 



qq and gg channels comprise about 300 and 800 one-loop 
diagrams, respectively. The most complicated ones are 
the 84 pentagons and 21 hexagons that contribute to the 
gg channel (see examples in Fig. [2]) and involve tensor in- 
tegrals up to rank five. Feynman diagrams are generated 
with two independent versions of FeynArts [33|, |3J], 
and one-loop amplitudes are reduced along the lines 
of Refs 



35l l36j using two in-house Mathematica pro- 



grams, one of which relies on FormCalc [37] for prelim- 
inary manipulations. The employed approach strongly 
mitigates the complexity inherent in Feynman diagrams 
by exploiting factorisation of colour matrices, reduction 
of helicity structures to compact spinor chains, and re- 
cycling of a multitude of common subexpressions. The 
treatment of rational terms of ultraviolet or infrared ori- 
gin is described in Appendix A of Ref. [35*1 . The reduced 
expressions are automatically converted into Fortran?? 
programs that evaluate colour/helicity summed quanti- 
ties with very high CPU efficiency. Tensor integrals are 
related to scalar integrals by means of numerical algo- 
rithms that avoid instabilities from inverse Gram deter- 
minants and other spurious singularities 31, 3^ . 



The presence of intermediate unstable top quarks in 
pp W+W~bb -|- X represents a non-trivial new as- 
pect as compared to previous NLO QCD studies of 
multi-particle processes. To regularise intermediate top- 
quark resonances in a gauge-invariant way we employ 
the complex-mass scheme [30| . In this approach the top- 
quark width Ft is incorporated into the definition of the 



renormalised (squared) top-quark mass, iif 



'irritTt 



In the on-shell scheme this complex parameter /i^ is iden- 
tified with the position of the pole of the top-quark prop- 



the tog-quark self-energy at p"^ 
Ref. m) 



fi^ via (see (4.25) in 



6^^, = ^ [S'-^(^?) + E*^i^(^?) + 2I]'.S(M?)] . (1) 



We note that an expansion of the occurring self-energies 



around the real point 



mt 



(as e.g. suggested in (4.2?) 



in Ref. [3C|) is not sufficient for NLO accuracy, because 
the top-quark self-energy is not analytic at the complex 
pole, p1 = Ht- The evaluation of one- loop scalar box in- 
tegrals in presence of complex masses represents another 
non-trivial aspect of the complex-mass scheme. In our 
calculation we employ the results of Ref. 38| , where ex- 
plicit analytic continuations have been presented for all 
kinematic box configurations that are relevant for physi- 
cal processes. 

The real corrections receive contributions from the 
2 — >• 5 partonic processes gg — >■ W+W~bbg and qq — >■ 
W+W~bbg, as well as from crossing- related gq and gq 
channels. The 2 — )■ 5 matrix elements are evaluated with 
Mad GRAPH [3^ and, alternatively, using the Weyl-van- 
der-Waerden formalism of Ref. [40|. To isolate infrared 
divergences and cancel them analytically we employ in- 
house implementations of the dipole subtraction formal- 
ism (4l| . Specifically this is done in dimensional reg- 
ularization with strictly massless light quarks (includ- 
ing b quarks) and alternatively in a hybrid scheme with 
small quark masses with the respective dipole subtraction 
terms from Ref. [42] 

Colour and helicity correlations that enter the sub- 
traction procedure are generated by means of Au- 
TODiPOLE [43| and, alternatively, in analytic form. To 



achieve sufficient numerical stability we perform the 11- 
dimensional phase-space integration using multi-channel 
Monte Carlo techniques with weight optimisation j3]. 
The integration of the dipole-subtracted 2 — >■ 5 contribu- 
tions is optimised by means of additional channels corre- 
sponding to the dipole kinematics. 

In the following we present predictions for the Teva- 
tron (pp collisions at 1.96 TeV) and the LHC (pp 
coUisions at ? TeV). In NLO (LO) QCD we employ 
MSTW2008NLO (LO) parton distributions ^ and 
describe the running of the strong coupling constant 
as with two-loop (one-loop) accuracy, including five 
active fiavours. Contributions induced by the strongly 
suppressed bottom-quark density are neglected. For the 
gauge-boson and top-quark masses we use rrit = 1?2 GeV, 
Mw = 80.399 GeV, and Mz = 91.1876 GeV. The masses 
of all other quarks, including b quarks, are neglected. In 
view of the negligibly small Higgs-mass dependence we 
adopt the Mh — ^ 00 limit, i.e. we omit diagrams involv- 
ing Higgs bosons. The electroweak couplings are derived 



from the Fermi constant = 1.16637 x 10' 



-5 GeV"^ in 



the G^-scheme, where the sine of the mixing angle and 
the electromagnetic coupling read = 1 



I and 
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a = ^/2G^M^sl,/n. For consistency, we perform the 
LO and NLO calculations using the top-quark widths 
Tt^Lo = 1.4655 GeV and Tt^NLO = 1.3376 GeV 0, 
respectively. Since the leptonic W-boson decay does 
not receive NLO QCD corrections we employ the 
NLO W-boson width Fw = 2.0997 GeV everywhere. 
Final-state quarks and gluons with pseudo-rapidity 
1 77 1 < 5 are converted into infrared-safe jets using the 
anti-fcx algorithm j^^]. For the Tevatron (LHC) we 
set the jet-algorithm parameter R — 0.4 (0.5) and 
apply the transverse-momentum and pseudo-rapidity 
cuts pT,b-jct > 20 (30) GeV, |?7b-jot| < 2.5. More- 
over, we require a missing transverse momentum 
of PT,miss > 25 (20) GeV and charged leptons with 
PT,i > 20GeV and |r/i| < 2.5. 

The LO and NLO W+W"bb cross sections at the 
Tevatron and at the LHG are plotted in Fig. |3] as a 
function of the renormalisation and factorisation scales, 
Mren = /^fact = A*. At the Tcvatron, where the cross sec- 
tion is dominated by the qq channel, at = mt we obtain 
(jT^v ^ 44.3i+J9-68fb and ct^cVq = 41.75;°;™ fb, where 
the uncertainties describe missing higher-order correc- 
tions estimated via scale variations in the range mt/2 < 
/i < 2mt. For the LHC, where the gg channel dominates, 
we obtain cr{;gc ^ 662.4t?^|:^ fb and a^^g = 840^7^ fb. 
Normalising the results to LO predictions at /x = m-t we 
obtain the relative NLO corrections K'^"" = 0.9421°;°°° 
and K^^'^ = 1.27tJ^;'j'i. The NLO corrections induce a 
moderate shift of the integrated cross section and reduce 
its scale uncertainty from about 44% (40%) to 9% (9%) 
at the Tevatron (LHC). This confirms the good conver- 
gence of perturbative predictions at the scale fj, — irit, a 
feature that is reflected also in the stable shape of the 
NLO curves in Fig. [3] 

To assess the impact of finite- width effects on the inte- 
grated cross section we have extrapolated our numerical 
results to the narrow- width limit Ft — ?> 0. In this region 
we observe a linear Ft-dependence, consistent with the 
cancellation of logarithmic soft-gluon singularities. At 
the Tevatron we find that finite-width effects shift the 
LO(NLO) cross section by about -0.8%(-0.9%). At the 
LHC we observe a qualitatively different behaviour: the 
shift induced by finite-width contributions is smaller in 
size and positive. At LO it amounts to -1-0.4%, and at 
NLO it becomes as small as the Monte Carlo statistical 
error (-^0.2%). 

To illustrate NLO and finite- width corrections to 
differential observables, in Fig. |4] we plot the invariant- 
mass distribution of a positron and a b-jet — the visible 
products of a top-quark decay — at the Tevatron. In 
narrow-width and LO approximation this kinematic 
quantity is characterised by a sharp upper bound, 
-^o+b — "^t ~ -^W' which renders it very sensitive to the 
top-quark mass. The value of rrit can be extracted with 
high precision using, for instance, the invariant-mass 
distribution of a positron and a J/^ from a i3-meson 
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cross sections at the Tevatron and the LHC. 



decay an observable that is closely related to M^+b. 
In Fig. |4]we clearly see, already in LO, small but non- 
negligible off-shell contributions that elude the kinematic 
bound. At NLO this feature becomes more pronounced, 
also due to QCD radiation that enters the b-jet without 
being emitted from its parent b quark. Below the 
kinematic bound we find very significant NLO effects. 
In the region 50 GeV < Mc+h < 150 GeV the shape 
of Me+b is strongly distorted, with corrections ranging 
from -1-15% to —30% (see lower plot). In the vicinity 
of the kinematic bound the NLO prediction is barely 
consistent with the LO uncertainty band. This example 
demonstrates the importance of 2 — > 4 NLO predictions 
for a precise description of the kinematic details of 
the W"'"W~bb final state and, more generally, for the 
top- physics programme at the Tevatron and at the LHC. 
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da r fb 1 pp -)■ W+(-i-!^ee+)W-(-j-/t-;y,,)bb+X 
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FIG. 4: Invariant mass Mg+j, of the positron-b-jet system at 
the Tevatron: absolute LO and NLO predictions (upper plot) 
and relative corrections w.r.t. LO at /i = mt (lower plot). The 
uncertainty bands describe mt/2 < fi < 2mt variations. 
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Note added: Shortly after the submission of this pa- 
per results of a similar calculation by the HELAC-OFF 
Collaboration have been published in Ref. (49| . 
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